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Abstract 

The stability of the flow of a viscous incompressible fluid over a curved compliant wall to 
longitudinal Gortler vortices is investigated. The compliant wall is modeled by a particularly 
simple equation relating the induced wall displacement to the pressure in the overlying fluid. 
Attention is restricted to the large Gortler number regime; this regime being appropriate to 
the most unstable Gortler mode. The effect of wall compliance on this most unstable mode 
is investigated. 
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I. INTRODUCTION 


Our concern is with the effect of wall compliance upon the stability characteristics of 
the flow of an incompressible fluid over a wall of variable curvature. Recent theoretical 
work on the effects of wall compliance on the stability of fluid flows has been restricted to 
flows over flat compliant walls (here by flat we are referring to the unperturbed state of 
the wall). This work has incorporated various models of the resulting wall displacements 
eg. viscoelastic walls 1 and spring backed walls 2,3 (both isotropic and anisotropic). Conse- 
quently, this work has considered various aspects of the stability of Tollmien-Schlichting 
waves and the effect that the underlying wall motion has on their stability characteristics. 
A review of these results can be found in Carpenter 4 . 

The effect of wall compliance on flows over curved walls has recently been considered 
experimentally by Yurchenko, Babenko & Koslov 5 . They restricted their attention to 
centrifugal instabilities, namely Gortler vortices. Although no precise measurements of the 
growth rates of the initial disturbances have been presented by these authors their results 
for a compliant wall, in which they plot regions of the Gortler number /wavenumber space 
corresponding to maximum amplification rates of the initial disturbances, show a shift to 
the left as compared to the results obtainable for the a rigid wall. These results then 
suggest that the effect of the compliant wall has a stabilizing effect on the growth rate of 
the longitudinal vortex structure. 

In the case of a rigid wall it has been shown by Denier, Hall & Seddougui 6 that 
there is a most unstable Gortler mode (ie a mode with maximum growth rate) in the 
parameter regime G 1 (here G is the Gortler number) and has wavenumber given by 
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k = 0(G 1/5 ) with growth rate 0(G 3/5 ). Such a result was found by considering both the 
inviscid Gortler modes, found in the parameter regime G 1 with wavenumber k = 0(1), 
and the right hand branch Gortler modes, found in the parameter regime G 1 and 
wavenumber k = 0(G ^ 4 ). By their very nature the inviscid modes do not feel the effect 
of the wall and so we expect that these modes are virtually unaltered by the presence of 
a compliant wall. In the case of the right hand branch modes it is known from the work 
of Hall 7 that these modes are confined to a thin region internal to the flow and away from 
the wall. Hence, we expect that these modes are also unaffected by the presence of the 
compliant wall. However, in the case of the most unstable Gortler mode it is known 6 that 
these modes are confined to a thin wall layer of thickness 0(G _1/5 ). We anticipate that 
these modes will be strongly affected by the presence of the compliant wall. For these 
reasons we will confine our attention to the most unstable Gortler modes. 

The outline of the paper is as follows. In §2 we formulate the problem under consid- 
eration and derive the governing equations for the fluid motion. The equation of motion 
for the wall displacement is assumed to have a particularly simple form (thus facilitating 
the analysis to be presented in §3); further refinements to the governing equation for the 
wall displacement can be made to take into account the effect of fluid substrates, elasticity 
of the compliant wall etc 4 . Such issues will not be pursued in this paper. Subsequently 
in §3 we will restrict our attention to the high Gortler number regime and consider the 
wavenumber regime k = 0(G 1 ^ 5 ). The eigenvalue problem for the most unstable Gortler 
mode will be derived and in §4 we present some results obtained from a numerical solution 
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of this problem. The implications of these results to the stability of flow over a curved 
compliant wall will be discussed. 


II. FORMULATION OF THE PROBLEM 

We consider the steady flow of a viscous incompressible fluid over a compliant wall 
of variable curvature. Assuming that L is a typical lengthscale over which the curvature 
changes, U ^ is a typical flow velocity a great distance from the wall, and v is the kinematic 
viscosity, we define a Reynolds number, Re, by 


Re = 


U X L 

v ’ 


( 1 ) 


and suppose that with respect to Cartesian axes x* , y * , z* the unperturbed position of 
the wall is defined by 

t r =LRe- 1 ' 2 g(x*/L). (2) 

Upon the imposition of some (suitably) small pressure perturbation the position of the 
wall is now defined by 

y* =LRe- 1/2 {y(x7L) + ^r?*(xVL,Re 1/2 z7L)}, (3) 


where 6r}* is the wall displacement whose motion is assumed to be governed by* 


LRe- 1/2 dBoVy = -6p\. 


( 4 ) 


Here Bo is the flexural rigidity of the compliant surface (we have assumed that the com- 
pliant surface is isotropic) and V 4 is the usual bi-harmonic operator 

v 4 = a 4 + 2d 2 , a 2 +a 4 .. 
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The term 6p\ is the small pressure perturbation. 


Defining nondimensional variables 

(x,y,z) = (x* ,Re 1/2 y* ,Re 1,2 z*)/L, 

(u + ,v+,u; + ) = (u*,Re lf2 v\Re 1/2 w , )/U oa , 

p + = p*/(pUl), p* = p* Q + 6p\, 
we restrict attention to the limit Re —* oo and write 

(u+, w+) = (a, V, 0) + 6(u,v, w) + 0(6 2 ), 

p+ = p(x) + 6Re~^ 2 p(x,y,z) + 0(6 2 ), pf = Re~ l/2 p. 

Here the basic velocity (a, V) depends only on x and y whereas the velocity perturbations 
( u,v,w ) depend on all three dimensionless coordinates. Substitution of these expressions 
into the Navier-Stokes and continuity equations, taking the limit 6 — » 0 with Re held 
fixed, equating terms of order <5°, (5 1 and taking the further limit Re — > oo in the resulting 
equations for (u,V,0,p) and (u,v, w, p) gives, to leading order in powers of i2e -1 , 

U a + Vy = 0, UU X + VUy = —p t + Uyy, (5) 


and 


U x + Vy + W x = 0, 


”1“ f)Uy “I" UU X “I" Vfly — Uyy | U|| , 

UV X + VVy + UV X + VVy = ~Py + Uyy + U„ , 
UW X + VWy = —P t + Wyy + Wjj, 

together with the governing equation for the wall displacement 


B 


d 4 fi 

dz 4 


= -p, 


( 6 ) 


( 7 ) 
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where B = Re 2 Bo/(pU^L 3 ) is a nondimensional constant which we initially take to be 

0 ( 1 ). 


By writing down the Taylor series expansion for the no-slip conditions about y = g(x) 


we obtain the boundary conditions 

0 = V = 0 y = g(x), 

u = — U y fj(x,z) t v = — E y rj(x, z), w = 0 
In the limit y — ► oo we have the requirement 


( 8 ) 


y = 9{x). 


U — » u t (x), (u,v,w)—*0 y — ► oo. (9) 

where u,(x ) is the dimensionless free stream velocity; in this paper we will confine our 
attention to the Blasius profile so that u,(x) = 1. In order to highlight the effect of wall 
curvature we make use of the Prandtl transformation 


y-*y + g , v-+v + g'n, v^v + g'u, 

whilst all other variables remain unchanged. The governing equations then become 

+ Vy = 0, UU Z + VUy = —p X + Uyy , 

U = V = 0 y = 0, G — > 1 y — > oo, 

and 

Ux Uy | — 0, 

UU X + VUy + UU X + VUy = Uyy + U , , , 

UV X + VVy + U'Ox + VVy + Gx(x)UU = —py + Vyy + U, , , 

UW X + VWy = — Pj + Wyy + W , Z > 

BfjzziX = Pi 


( 10 ) 


( 11 ) 


u — —Uyfi, v = w = 0 y = 0, 
(u,v,w) -> 0 y — ► oo. 
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In (11) we have replaced 2 g xx by Gx(a;) where G and \ will be refered to as the Gortler 
number and wall curvature respectively. 


Finally taking the Fourier transform of the system (11) in z, eliminating w, p and the 


wall displacement rj from the resulting equations gives 


O 

Uu m + VUy + uTL m + vUy = Uyy — k Uy 

{ ^ + fc^Dy } V + d X Uyy + { tZ #a y + ft X + \Gu} U 

f _ d 2 l2 \ f d \ (12a,6) 

+ | Uyy U ^2 + k Q j V X + 2 | S®y + Uj, ^ J U* 

Vyyyy flVyyy { Uy 21/ } Uyy “|“ { Sjy “|“ fc Uy 0, 

together with the boundary conditions 


u — ^ (^*vv d* u vvy ) d" A; u® •+■ u{u xx -t* u®y) -h *0{u X y + Uyy)} y — 0, 


v = v = 0 


V = 0 , 


(13) 


u,v, v — > 0 y — > oo. 

Here A: is the transform variable, (u, u) are the transforms of (u, v). For 0(1) values of the 
wavenumber k and finite values of the Gortler number G the solution of the system (12,13) 
can only be found by numerical integration. However, in the asymptotic regime G 1 
progress can be made analytically and it this parameter regime which we consider in the 
following section. From the boundary condition (13) we see that the case of the rigid wall 
is recovered in the limit B — » oo. 


III. THE MOST UNSTABLE GORTLER MODE 

In the case of a rigid wall it has been shown 6 that the system (12,13) has a mode with 
maximum growth rate in the parameter regime G 1 and the wavenumber of this most 
unstable mode is k = 0(G 1 / 5 ). Such a result was found by considering the structure of 


6 



both the problem of the inviscid Gortler modes, which are found in the regime k = 0(1), 
G 1, and that of the right hand branch modes, found in the parameter regime G 1 
and fc = 0(G 1/4 ). 

By inspection of the structure of both the inviscid modes and the right hand branch 
modes 9 it is clearly seen that in the case now under consideration these modes are essen- 
tially unchanged from those relevant to the rigid wall problem. In fact, an analysis of the 
system (12,13) shows that in these two parameter regimes the effect of wall compliance 
is encountered as a lower order effect and does not contribute, at leading order, to the 
structure elucidated in Ref. 6. This result is easily seen if one considers the position of 
both the inviscid and right hand branch modes; the inviscid modes (by definition) do not 
feel the effect of the wall whereas the right hand branch modes are situated at a position 
internal to the flow regime and so will not feel the effect of the compliant wall. 

For these reasons we will restrict our attention to effect of wall compliance on the 
most unstable Gortler mode. To consider this mode, we note 6 that this mode is confined 
to a wall layer of thickness 0(fc -1 ) where the wavenumber scales as k = 0(G^ 5 ). We then 
write 

<p = ky, k = AG 1 / 5 . (14) 

We expand the mean flow in this layer as 

Q ~ n(x)k~^<p + • • • , 
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and expand the velocity field (u, v) and the surface displacement ij as 


/ * \ 

u={u 0 +G- 1 ' 5 U 1 + ---}exvyG 3 ' 5 J j3(x)dx j , 

( X \ 

v =G 2,b { V 0 + G~ 1/5 Vi + • • •} exp ^ G 3,b J f3(x)dx j , (15) 

rj = {770 + G~ 1,5 T]i H } exp yG 3,b J /3(x)dx j , 

where we have anticipated growth rates of 0(G 3,b ) (see Ref. 6). Substitution of (14) and 
(15) into the governing equations (12a, b) we obtain, to leading order, 

r d 2 nP<P t X ,, H T/ 

U °=T2 Vo ' 


d^ A 3 

wVo 


(16) 


d(p 2 “‘j - a 3 

An evaluation of the boundary conditions (13) shows that the effect of wall compliance 
first becomes important when the compliance parameter B = 0{G ~ 1|/5 ). Hence, we define 


B 0 = BG 1,b . 


The boundary conditions appropriate to (16) are then found to be 

u, d 3 Vn 1 

Po = 5^TO’ Ko = 0 ’ Fo = 0 ” = 0 ’ 

Uq = 0, Vo = 0, Vq = 0 <p = 00 . 

As noted above, the case of the rigid wall is recovered in the limit Bo — > 00 . 

For computational purposes it is convenient to eliminate the parameters /x, x from 
(16) by writing 



A = (XP 2 ) 1/5 A, £=(xV /5 /*, B 1 = Bo{xh 2 ) 1 '\ U 0 = pUo, Vo = (x^ 2 ) 2/5 Vb, 
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in which case (16,17) may be rewritten as 

,U-, _Vo 

l dip! A3 1 J t ' 0- A2’ 

j2 




<^0 


#0 = 


1 d 3 Vo 

A 3 Si dip 3 ’ 


(18) 


Vo = v-n = 0 


v? = 0, 


Uq = Vo = Vq = 0 <p = oo. 


For a given value of the parameter B\ the system (18) constitutes an eigenvalue problem 
for the growth rate (3 = /?(A). 

Before proceeding with a discussion of the results of our numerical integration of the 
system (18) we note that the limiting inviscid form of the governing equations is recovered 
by taking the limit A — ► 0 while the limiting form of the right hand branch governing 
equations are obtained by taking the limit A — > oo. In fact, from (18) it is easily seen that 


(3 ~ A 1/2 , A — > 0, /? A 2 , A — > oo, 


thus demonstrating that for some intermediate value of A there exists a most amplified 
mode. 


IV. RESULTS AND DISCUSSION 

The system (18) was solved numerically using a finite difference scheme in the in- 
dependent variable <p. For a given value of the parameter B\ and a fixed value of the 
wavenumber A an initial guess for the growth rate (3 was supplied and the program was 
run until convergence was achieved. The eigenvalue problem (18) was then solved by 
continuously iterating in the wavenumber A to generate eigenvalues (3 = /3(A; B\). This 
procedure was repeated for various values of the parameter B\. 


9 



For a given value of B\ the system (18) has an infinite number of eigenvalues /?„ = 
/3 n (A; B\). In Figure 1 we present a plot of the first eigenvalue (3 as a function of A for 
various values of the parameter B i; these eigenvalues correspond to the most unstable 
Gortler mode. In Table I the maximum eigenvalue /? ma.v is given together with the corre- 
sponding critical wavenumber A cr - We see that the introduction of wall compliance has a 
stabilizing effect in that the maximum value of the growth rate is reduced. From Figure 1 
we see that as the wall compliance is increased (ie B\ is decreased) the maximum growth 
rate decreases until a threshold value is reached. Below this threshold value of Bi the 
stabilization due to the wall compliance remains constant; from Table I we see that the 
stabilization is 5% compared to the case of the rigid wall. 

From Table I we see that as the compliance parameter is decreased from the rigid wall 
value there is a slight increase in the critical wavenumber until a threshold value of the 
parameter Bi is reached after which point the critical wavenumber decreases in magnitude. 

In figures 2 and 3 we present the eigenfunctions for the maximum eigenvalues presented 
in Table I. The effect of wall compliance is to decrease the amplitude of the eigenfunctions. 
Otherwise there is little qualitative change from the case of the rigid wall. 

In conclusion, the presence of a compliant wall has the effect of stabilizing the longi- 
tudinal vortex motion which is in qualitative agreement with the available experimental 
results 5 . However, this stabilization is relatively small and hence we may conclude that 
the stability of the flow over a curved compliant wall is essentially unchanged from that of 
the case of a rigid wall. It remains an open question as to whether either finite amplitude 
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effects or the presence of wall modes (ie modes induced by the motion of the compliant 
wall) will significantly influence the structure discussed in this paper. 
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Table I. Maximum growth rate and critical wavenumbers. 


B i 

/^max 

^cr 

oo°) 

0.3135 

0.475 

100 

0.3113 

0.49 

10 

0.3039 

0.513 

1 

0.2985 

0.487 

0.01 

0.2977 

0.476 


°)The case of the rigid wall. 
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